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Abstract 

Single-particle kinetic energy spectra and two-particle small angle correlations of 
protons (p), deuterons (d) and tritons (t) have been measured simultaneously in 
61A MeV ^^Ar -|- ^^Al, ^^^Sn and ^^^Sn collisions. Characteristics of the emission 
sources have been derived from a "source identification plot" (Psource-EcM plot), 
constructed from the single-particle invariant spectra, and compared to the com- 
plementary results from two-particle correlation functions. Furthermore, the source 
identification plot has been used to determine the conditions when the coalescence 
mechanism can be applied for composite particles. In our data, this is the case only 
for the Ar + Al reaction, where p, d and t are found to originate from a com- 
mon source of emission (from the overlap region between target and projectile). In 
this case, the coalescence model parameter, pQ - the radius of the complex particle 
emission source in momentum space, has been analyzed. 
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1 Introduction 



A large body of studies, both experimental [1-8] and theoretical [9-17], has demon- 
strated that intermediate energy heavy ion collisions proceed through a complicated 
reaction mechanism. For a wide range of impact parameters, ranging from periph- 
eral to near central, the mechanism is dominated by dissipative binary collisions, 
with early dynamical emission followed by statistical evaporation. In particular, at 
E/A ~60 MeV, the emission of light particles originates from (at least) three sources 
[2,18], a quasi-projectile source (QP), a quasi-target source (QT), and an intermedi- 
ate velocity source (IS). The IS accounts for dynamical emission, described by early 
nucleon-nucleon collisions and by other preequilibrium processes, such as emission 
from a low density neck region [18-30] . Since the particles emitted from these sources 
overlap in angle, energy and emission time, it is not possible, on a particle-by-particle 
basis, to identify which source each particle came from. 

A powerful experimental procedure to gain insight into the reaction mechanism is 
two-particle interferometry, used to determine the space-time extension of the parti- 
cle emitting sources [31,32], and to extract the emission time sequence of neutrons, 
protons and their composites [33] . Our recent interferometry analysis of Ref . [34] on 
the 61 MeV/nucleon ^^Ar -|- ^^Al reaction, has confirmed a dissipative binary reac- 
tion scenario with an important dynamical emission component from the IS created 
in the overlap zone. In particular, an angular dependence of the correlation func- 
tion has emerged. The stronger pp correlation observed at backward angles, indicates 
emission from a shorter lived source. This has been interpreted as an enhancement 
of the early dynamical emission component seen at backward angles. The average 
emission time of deuterons from this source has been found to fall in-between those 
of neutrons and protons, as expected if deuterons are formed mainly by coalescence 
[35]. 

To complete and complement the interferometry analysis, the information contained 
in the single-particle energy spectra is explored in this paper. In particular, we inves- 
tigate the structure of light particle emission sources as it follows from the kinematic 
features of the single-particle {p, d, t) energy spectra. A common approach to study 
single-particle energy spectra, is to fit the data to a moving source parameterization, 
using multiple sources to describe preequilibrium and evaporation sources. Contribu- 
tion from three sources is normally assumed: a QP, an IS and a QT [18,36-39]. While 
such a phenomenological approach can give a reasonably good description of double 
differential particle cross sections, the dynamics of the nuclear interaction and the 
mechanism of particle production remains hidden in the large number of fitted pa- 
rameters. This description is further complicated by the fact that the three-sources 
relative yields vary with impact parameter, so that the parameters utilized in the 
description only represent average values. In this paper we investigate an alterna- 
tive approach to the fitting procedure with three moving sources. This approach is 
feasible when a high-precision energy calibration is available, as is the case in the 
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present data set [40]. In this approach, we construct a Psource^Ecu plot from the 
invariant energy spectra. Here, f3source is the source velocity in units of c, and Eqm is 
the particle energy in the source frame. Inspection of this plot reveals a continuum 
of source velocities, in accordance with the expected trends for a dissipative binary 
collision mechanism of heavy ions. For the Ar + Al reaction, the continuum of source 
velocities is found to be common for p, d and t charged particles, and therefore we 
can proceed to apply a coalescence mechanism, which is proved to be the dominant 
process of dand Remission from the overlap region between target and projectile. The 
coalescence model parameter, pQ - the radius of complex particle emission source in 
momentum space, is analyzed, which provides valuable complementary information 
to the two-particle correlation size measurement. 

The paper is organized as follows. Sec. 2 presents the details of the experiment. Sec. 
3 presents the source analysis for the ^^Ar + ^^Al reaction. The angular dependence 
of p, d and t single-particle energy spectra is presented in Sec. 3.1; the comparison 
with results from Firestreak model calculations highlights some simple physics behind 
the emission mechanism. Sec. 3.2 introduces the jdsource^EcM plot, and discusses the 
wealth of information it contains on the multiple source reaction scenario. In Sec. 
4, the information extracted from the fisource^EcM plot is compared with the results 
from our previous interferometry analysis [34]. Coalescence analysis of the d and t 
spectra is performed and discussed in Sec. 5. Sec. 6 extends the Psource^EcM source 
analysis to 61A MeV ''^^Ar + ^^^'^^^Sn coHisions. A summary can be found in Sec. 7. 
Appendix A describes the details of the Firestreak model calculations. Appendix B 
describes the details of the construction and interpretation of the Psource^EcM plot. 



2 Experimental details 

The experiment was performed at the AGOR cyclotron of KVI (Groningen). The E/A 
= 61 MeV 36^j^i4+ p^iggcl beam (~ 0.5 -^1.0 nA) impinged on targets of ^''Al, ^^^Sn 
and ^^^Sn (1.8 mg/cm^ thick). The beam intensity was monitored by a Faraday cup. 
The experiment measured simultaneously single-particle energy spectra and small 
angle two-particle correlations, in coincidence with forward emitted fragments. The 
requirement of at least one fragment detected in the forward direction biases our 
collected data towards semi-peripheral coUisions [41]. 

The particles were detected using three complementary multidetector systems: the 
EMRIC charged particle detector array [40,42], the EDEN neutron detector array 
[43,44], and the KVI Forward Wall (FW) phoswich detector [45]. The experimental 
data were collected in two modes. In the first mode, both EMRIC and EDEN reg- 
istered inclusive particle spectra in coincidence with the FW. In the second mode, 
coincidence events between EMRIC, EDEN and EMRIC*EDEN elements were col- 
lected in coincidence with the FW. 
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The EMRIC charged particle detector array consisted of 16 CsI(Tl) elements, each 
one subtending a solid angle of 3.6 or 5.1 msr. The detector modules were arranged 
in two groups, consisting of single elements and triple elements. The complete ge- 
ometry, optimized for the interferometry requirements, is presented in Ref. [44]. All 
detectors were operated in air, outside a scattering chamber of (3 mm thick) stain- 
less steel and 28 cm in radius. 100 fim thick Capton foil windows were installed in 
front of the EMRIC elements. One extra element, equipped with a ^^^Am a-source 
and placed at 30° was used to monitor the shift in the a-peak position caused by 
the variation of beam intensity and, consequently, the 7-background, detected by all 
EMRIC elements. Temperature monitoring of the CsI(Tl) elements was also applied 
[40]. The detectors and measuring technique are described in details in Ref. [40]. The 
i,2,3jj isotopes were identified via the pulse-shape analysis method. An accuracy < 1% 
was achieved in the energy calibration for each of the 100 mm long crystals by using 
the AE(Si)~E(CsI(Tl))/PMT method. The validity of the calibration was checked 
by comparing the particle spectra measured by detectors located at the same polar 
angle. Isotopic resolution was obtained for Z—1 starting from energy threshold ~ 8, 
11 and 14 MeV, for p, d and respectively. 

Two important corrections, both energy dependent, were introduced in the measured 
particle spectra. These corrections are quite common for experiments with long scin- 
tillators. First, the correction on the loss of charged particles due to their inelastic 
interaction with the scintillator material. Experimental data and calculations [46] 
may be approximated with a good precision by the equation, 

hnei = 1 + h{Z, A)E[l - exp{-E/im% (1) 

where E is the particle kinetic energy in MeV, A) is a fitting parameter. 

Second, the correction on the change in real solid angle, ^reaii of the scintillator along 
its length, 

ko, — ^input/^real = (H — ^-E^'^^)"^. (2) 

Here Qinput is the solid angle defined by the front face active area of the scintillator, 
L(cm) is the distance from the target to the scintillator front face, and c(Z, A) is a 
fitting coefficient, extracted from the energy-range relation for particle (Z, ^4) in Csl 
crystal [47] ^ . The measured spectra have been corrected by the factor k = kind ■ k^, 
i.e. N {actual, E) = k ■ N {measured, E). As an example, the total correction factor 
for protons is calculated to 1.76 for E=200 MeV. 



^ The coefficients b{Z, A) and c{Z, A) are equal to 0.0014, 0.00090, 0.00066 and 0.00140, 
0.00155, 0.00115 for p, d, t, respectively. 



4 



3 Source analysis for E/ A = 61 MeV ^^Ar + ^^Al collisions 



3.1 Single-particle energy spectra 

Proton, deuteron and triton energy spectra, produced in 61 MeV/nucleon ^^Ar + ^^Al 
collisions, were measured in the energy intervals ~ 8(ll)(14)-250 MeV for 8 angles 
in the polar angle region 30°-114°. Fig. 1 presents proton (upper panel), deuteron 
(middle panel) and triton (bottom panel) data for some laboratory angles, as listed 
in each panel. 

We have compared the proton data (Fig. 1, upper panel) to the results of a Fire- 
streak model calculation (histograms) and normalized to the data at 78°. The ba- 
sic assumptions of the Firestreak model are reviewed in Appendix A. This model 
uses a quite drastic approximation of the geometry and treats the overlap region of 
the target and projectile as created in a totally inelastic process. Yet, a reasonably 
good description of the angular and energy distributions of the protons is achieved. 
It should be stressed that only one parameter - the maximum source temperature, 
7"source= 13.2 MeV is introduced. The quite good agreement with the experimental 
data shows that the nature of the binary dissipative collisions is closely related to the 
participant-spectator picture of the interaction (see also Sec. 3.2). 

3.2 Source identification, Psource-EcM plot 

The simplest way to describe the particle spectra is to adopt the participant-spectator 
description and assume emission from three thermal moving sources. For each emit- 
ting source, the velocity, (^source, source temperature, Tgource, Coulomb barrier, V^, 
and one normahzation constant are treated as free parameters. At least 12 parame- 
ters must be determined by a simultaneous fit to the particle spectra at the measured 
angles. To overcome the problem of coupling between the parameters, some of them 
must be kept fixed. Such a parametrization gives a reasonably good description of 
the particle spectra and angular distributions [18,36-39]. 

As an alternative way we construct a "source identification plot" from invariant cross 
sections. We obtain the source velocity, ^source, and the kinetic energy of the particles 
in the source frame, Ecm, from invariant cross sections, using different lab energies 
and different angular combinations. In this way, we construct a Psource^EcM plot as 
if only one source was responsible for emission of particles. When several sources are 
present, the extracted value of ^source represents an effective average source. However, 
and this is an important point, the effective source will have different characteristics, 
depending on the angles of the chosen detector pairs. By careful analysis of the 
source identification plot, the properties of the contributing sources may be revealed 
to a certain extent. In Appendix B we discuss the derivation and properties of such 
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Psource^EcM plots further. Fig. 2 presents the results from the Psowce^EcM method 
with data from ^^Ar + ^^Al colhsions. Protons (full symbols), deuterons (big open 
symbols) and tritons (small open symbols) have been studied. The results are shown 
for some representative pairs of angles (selected for their highest statistics, taken with 
the triple EMRIC detectors). The value of Ecm is plotted in MeV/nucleon. The sohd 
lines are drawn to guide the eye. Values of the source velocity corresponding to total 
fusion {^source = 0.205) and to the projectile velocity {jSsource — 0.361) are marked by 
the arrows on the y-axis. 

This plot contains various types of information, which we now proceed to discuss. 

1) A main feature is that, in the intermediate region of Ecm, there is a clear spreading 
of values of Psource- This is due to the contribution from different sources (QT, QP 
and IS) and to the broad impact parameter range of our apparatus, leading to a 
continuum of source velocities, covering the region 0.07 < (^source <0.22 (common to 
all charged particles). This is in accordance with the expected trends for a dissipative 
binary collision of heavy ions. The observed continuum of emission sources is quite 
naturally reproduced in the nuclear Firestreak model description discussed in Sec. 
3.1 and in Appendix A (Fig. 4). 

2) One can identify a limiting region at high Eqm, where the source velocity converges 
towards the value l3source ^ 0.21. This indicates that the particles emitted at different 
angles, with the highest kinetic energy, originate from a source (IS) with a well defined 
velocity. 

3) Another limiting region can be identified at low Ecm and backward angles. Here 
the different lines tend to converge towards the QT velocity. These low kinetic energy 
particles may be identified as emitted by the QT source (compare with Fig. 5) . 

4) The "bump" observed for low Ecm and forward angles (the 30°/54° pair), repre- 
sents a third limiting region. This "bump" indicates that particles with a quite low 
kinetic energy are emitted by a source with a high velocity. These particles may be 
identified as emitted by the QP source (compare with Fig. 5). This region is not 
completely explored by our experimental apparatus, due to the limited coverage at 
forward angles. Notice that the QP "bump" disappears if the source identification 
plot is constructed for particles detected in coincidence with 4 or more fragments in 
the FW, a condition that biases the collected events towards less peripheral collisions. 

5) The J), d, t particles are emitted by common sources (i.e. the corresponding points 
lie on the same line in the plot). This is a key point that justifies the application of 
the coalescence mechanism of d and t production (see Sec. 5) and, maybe, a unique 
lifetime for their emission. 
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4 Comparison with interferometry results for E/A = 61 MeV ^^Ar + 
^^Al collisions 

The source characterization from the ^source — Ecm study, can be compared with the 
interferometry analysis carried out in Ref . [34] . There it was found that the correlation 
strength of pp pairs measured in the "forward" angular region (54°/30°) is smaller 
than the correlation strength measured backwards. This indicates a larger space-time 
dimension of the effective source that emits protons in the forward direction. This is in 
agreement with the Psource — Ecm result that detectors at "forward" angles, 54°/30°, 
register particles originating from a superposition of two sources, the IS and the QP 
source [11]. The QP decay reveals itself as the "bump" at Eqm ~ 25 MeV/nucleon. 

In contrast, an effective source of smaller space-time size was deduced from the 
stronger pp correlation function measured at backward angles [34] . This result seems 
to be in contrast with the source identification plot in Fig. 2, where at backward angles 
(78°/102°) a QT can be identified at low Eqm energies. However, one should notice 
that in the correlation function, the yield from the different sources (an information 
not seen in the ^source — Ecm plot) plays an important role. Since the measured yield 
from the QT source is low, the dominating contribution to the correlation function 
actually comes from the region of high Ecm energies in the source identification plot, 
dominated by the IS. 

Finally, also the results of the pp correlation functions with a high total momen- 
tum cut {Ptot > 450 MeV/c) are in agreement with the information derived from 
the Psource^EcM plot. In Ref. [34] it was found that the high total momentum gated 
pp correlation functions are enhanced and similar in size at backward and forward 
angles, indicating emission from an IS of small space-time dimension. The pp correla- 
tion function measured forward is more substantially enhanced. This is because the 
momentum gate at forward angles effectively corresponds to suppressing the contri- 
bution from the QP source, leaving the IS at high Ecm (high Ptot) as the dominating 
source. At backward angles, on the other hand, the effect of the high total momen- 
tum gate is only a small enhancement, since the contribution from the QT source 
was already small (see discussion above). 



5 Coalescence analysis for E/A = 61 MeV ^^Ar -|- ^^Al collisions 

The common sources for d and t emission. Fig. 2, suggest their common origin, f.e., 
through the coalescence mechanism. Calescence models are reviewed in Ref. [35] . 

An empirical nucleon coalescence model assumes a simultaneous emission of protons 
and neutrons (whose relative momentum is less than a certain value po) by the source 
to form a complex particle {Z,A). Indeed, our study of the particle emission time 
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sequence via two-particle correlation functions [34] is consistent with this assumption 
(as discussed in Sec. 1). The invariant momentum spectrum (or momentum space 
density) of the composite particle can be expressed through the neutron and proton 
invariant spectra, 

EAd^UA/d^PA = BAiEpd'ap/d'ppf ■ {End^ajd'pnf ■ (3) 



Here p^=App^„. Ep, En and Ea are the proton, neutron and composite particle total 
energies ^ . The "coalescence factor" is given in the classical approximation as: 



Here m is the nucleon mass, Sa is the spin of the particle, (Tq is the "hot zone" 
formation cross section, and po is the "coalescence radius" in momentum space. The 
value of (7o cannot be defined unambiguously. Usually, it is taken as the geometrical 
reaction cross section. The choice of ctq influences the extracted po value. 

If chemical and thermal equilibrium are achieved in the "hot zone" volume [48], po can 
be related to the volume Vq within which the nucleons coalesce to form a composite 
particle, 

V, ~ (5) 



For the dissipative binary type of collisions studied here, charged particles are emitted 
mostly in preequilibrium processes (with emission time < 50 fm/c [11,34]). Thus we 
take pq only as a parameter of the coalescence model. 

Coalescence model calculations of d and t spectra are presented in Fig. 1 for some 
representative angles. The histograms drawn for d (middle panel) and t (bottom 
panel) represent the calculation in coalescence model with po= 100 and 135 MeV/c 
for d and respectively. Agreement between the model calculation and the data is 
very good, with the exception of the 30° angle, where the low-energy part of the 
spectrum comes from the QP decay. (It should be noted here that the coalescence 
analysis must be restricted to the kinematical region for which particle emission from 
the QP and QT sources is negligible). 



^ Because of experimental restrictions, the spectra of protons and neutrons were registered 
by our apparatus at slightly different angles. Therefore, in the coalescence analysis we 
assume that they are equivalent for a given angle, but scaled by a factor Rn,p = {Nt + 
Np)/[Zt + Zp), corresponding to the ratio of neutron to proton numbers in the projectile 
and target. 
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Quantum mechanical treatment of the coalescence mechanism [49] assumes that the 
momentum distributions of the emitted particles can be described by density matri- 
ces. The coalescence volume is related to the internal wave function of the composite 
particle and of the spatial (Gaussian) distribution of nucleons in the source. Notice 
that the model does not impose thermal or chemical equilibrium. The coalescence 
radius is expressed as: 

p,{A)^const{N,Z,A,S)[^ + ^Y/', (6) 

where Ra and Rq are the radii of the composite particle and of the emitting source, 
respectively. The source radius extracted by applying the density matrix coalescence 
model is ~ 3.5 fm, both for deuterons and tritons. 



6 Source analysis for E/A = 61 MeV ^^Ar + ^^^'^^^Sn collisions 

Fig. 3 presents the source identification plot Psource-EcM constructed from the in- 
variant spectra of protons, deuterons and tritons emitted in 61A MeV ^^Ar -|- ^-"^^Sn 
(upper panel) and ^^Ar -|- ^^^Sn (lower panel) coUisions. 

Inspection of the plots reveals that: 

1) A source velocity continuosly varying with Eqm is again extracted (as for the Ar 
-I- Al system), in the region of 0.03 < Psource <0.2. 

2) The results for the two Sn targets are very similar, which indicates very little 
isospin dependence in the source velocities. A somewhat larger contribution of the 
target residue is noticeable at backward angles in the Ar -|- ^^^Sn system. 

3) The Ar + Sn results do not show a QP limiting region. As compared to the Ar + 
Al plot, the "bump" observed in the Psource^EcM dependence for low Ecm, forward 
angles, is absent in Ar -|- Sn. This indicates that our experimental apparatus misses 
the QP component in the direct kinematics Ar -|- Sn reactions, due to lack of forward 
angle coverage. 

4) The sources are common for p, d and t at high energies {Ecm >30 MeV/nucleon). 
In the low Ecm range, instead, protons deviate from deuterons and tritons. Inspection 
of this deviation reveals that the QT component is more sizable for proton emission. 
This implies that a coalescence model analysis of the deuteron and triton production 
is not feasible in Ar + Sn, since, as it was pointed out in Sec. 5, the coalescence 
analysis must be restricted to the kinematical region for which particle emission from 
the QP and QT sources is negligible. 
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7 Summciry 



We have proposed a novel and simple method to analyze single-particle energy spec- 
tra. This method allows to characterize the multiple sources of particle emission that 
are present in dissipative binary, intermediate energy heavy ion reactions. The method 
consists in constructing a source identification plot f3 source- Ecm- If the accuracy of the 
energy calibration is sufficiently high, this plot can be generated using the invariant 
cross section measured at different laboratory angles. 

Inspection of this plot suggests that 61 MeV/nucleon Ar + Al, Sn collisions are 
characterized by sources with continuously varying velocities, that tend to converge 
towards a well defined source velocity for particles emitted with high Ecm- Thus, 
this analysis supports a dissipative binary reaction scenario, with an important con- 
tribution from an intermediate velocity source. 

The results from the Ar + Al collisions have been compared with the two-particle cor- 
relation function analysis performed for coincidence data from the same experiment 
[34] . A connection between the apparent source size (containing space-time contribu- 
tions from different sources) extracted from the correlation functions and the results 
from the Psource-EcM plot has been established. 

The Psmirce-EcM plot also allows to compare the origin of emission of particles of dif- 
ferent types, and to reliably establish whether a common origin, necessary condition 
for the applicability of a coalescence mechanism, is met. In our data, we have found 
that this is the case for the Ar + Al collisions, but not for Ar + Sn. Application 
of the coalescence model to the Ar -|- Al data, has allowed to extract a coalescence 
radius, common for deuterons and tritons. 
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APPENDIX A : Firestreak Model 

The Firestreak model introduces a possibility to describe the fluctuations in souce 
size and excitation energy of the strongly interacting fireball. In our approach we 
follow essentially the prescription given in Ref . [50] . A nuclear density distribution is 
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incorporated for both colliding nuclei. The collision is treated as a totally inelastic 
process in the overlap region of the target and projectile. The interaction proceeds via 
coUinear streaks of nuclear matter of both nuclei. Each of the streaks is characterized 
by the value of "projectile fraction", 

rj ^np/{np + nt), (7) 



where np{nt) is the number of contributing nucleons from the projectile (target). The 
local value of source velocity, 

P{V)^V- Pbeam: (8) 



and the internal energy per nucleon, 

t{ri) = ri{l - ri)tbeam, (9) 

are both expressed through the value of rj. Here Pbeam is the velocity of the projectile, 
and theam is the kinetic energy per nucleon of the projectile. The geometrical aspects 
of the interaction (in units of cross section) are contained in the "yield function", 

Y{n). 

Fig. 4 presents the results of our numerical calculation for ^^Ar + ^''Al collisions at 
61A MeV. We used a three-parameter Fermi-type nuclear density distribution p(r) for 
both colliding nuclei. The parameters are taken as compiled in Ref. [51]. A parabohc 
interpolation was used to get parameters for the nuclei of interest in this work. The 
bins of Y{r]) near r] ^ correspond to contributions from the target spectator, with 
quite low internal (excitation) energy and source velocity. The bins of Y{r]) near f] ?a 
1 correspond to contributions from the projectile spectator, with quite low internal 
energy but with high source velocity. 

The Lab momentum space density for particle of type j is given by the expression, 

n 

FAp) =J:YiV^K'^M^^)] ■ fAp\tiV^)], (10) 
1=1 

where Jp'^p is the Jakobian transformation from the CM to Lab frame, and fj[p\ t{r]i)] 
is the CM frame momentum distribution for protons. Since we apply here only a phe- 
nomenological description of the proton spectra, the choice of fproton[p :t{Vi)] is quite 
arbitrary. Furthermore, we do need the assumption that particles in the local source 
have come to thermal equilibrium in their CM frame. So we accept a classical Maxwell 
momentum distribution in the CM frame of the emitting system, 

fproton{p,t) ~ {l/t^)exp{-^^^), (11) 
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with a surface proton emission. The Coulomb barrier is taken in the source frame, 
Vc=2.0 MeV. Results from the calculations are presented in Fig. 1 (protons). To get 
consistency with data, the "effective" temperature, 2t{r]i)/3, is increased by factor 
a— 1.32. If one postulates, as in the hydrodynamical approach [52], that the nuclear 
fluid is viscous, then friction arises over the region of overlap of the target-projectile. 
Work done against this friction force is transformend into a thermal energy, which 
essentially increases the participant and spectators excitation energy [50] and this 
motivates the a-factor. 

APPENDIX B: Source velocity reconstruction by the Psource~EcM niethod 

To facilitate the reconstruction of the sources we use different detector combinations 
to extract source velocity, ^source, and particle energy in the source frame, Ecm, from 
the invariant spectrum. If there is only one unique emission source from a reaction, 
the invariant cross sections at different angles can be used to determine the source 
velocity. The Lab energies Ei and E2 of particles emitted at different angles with equal 
invariant cross sections should thus coincide in energy (Ecm) after a transformation 
to the source system, moving with ^source in the laboratory. This allows to extract 
the value of ^source and Ecm from the shift in particle kinetic energy, Ei — E2, for 
each pair of Lab angles 6*1 and 6*2 using two equations, 

Ecm — El — 2(Ei ■ EgourceY^'^ ■ COsOi + Egource- (12) 
Ecm = E2 — 2{E2 ■ EgourceY^"^ ■ COSQ2 + E source- (13) 

Here Esource=i^l3lource/'^ the Lab kinetic energy of particle m, having the velocity 
^source- In fact, the described procedure means extracting Ecm and ^source by restoring 
the hypothetical invariant cross section in the source frame by using Eqs. (12, 13). It 
should be noticed that for this method to yield reliable results, the energy calibration 
must be known with a high accuracy. For example, an accuracy ±3% in the energy 
{Ei,E2) yields approximately ±10% error in 

By combining detectors at different laboratory angles and using particles emitted with 
different energies a number of different combinations of (^source-EcM can be obtained. 
The method described above assumes only one source present. When several sources 
are present, the extracted values of P source and Ecm represent an effective or average 
source. 

To correctly interprete the ^source and Ecm obtained for our experimental data, we 
construct a ^source and Ecm from a known input. For this purpose we use the three 
source fit by G. Lanzano et al. [18], which has been obtained from the double dif- 
ferential cross sections for protons emitted in E/A = 60 MeV '^^Ar + ^^Al collisions 
measured in the angular range 12.5°-112.5°. A standard way to describe the proton 
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spectra and their angular dependence has been apphed, by assuming particle emission 
from three thermal moving sources. In the source reference frame, particle evapora- 
tion is described as a volume emission of Maxwell-Boltzmann type. For each source, 
the values of Psource, "^source (MeV), Vc (MeV) and a normalization constant, a (in 
barn), are treated as free parameters. The corresponding fit parameters for 3 types 
of moving sources, QT, IS (denoted QF in [18] and in Fig. 5) and QP are determined 
by a simultaneous fit to the proton spectra at the 13 measured angles. The fit 
parameters, as deduced by the authors [18], are listed in Fig. 5. 

We used these parameters to calculate the invariant proton spectra. We then applied 
the procedure described to reconstruct ^source- The results are presented in Fig. 5. 
The upper panel presents the results when particle emission from one single source is 
assumed (QP, IS and QT are illustrated in the same plot). The figure shows that the 
reconstructed (3 source is in general equal to the input values, and the small deviations of 
some points illustrates the accuracay of the method. Another important obeservation 
is the range of energies populated by the different sources. Because of kinematical 
reasons the fast moving QP cannot populate low energies ^ . This is important to 
have in mind to correctly interprete the "effective averaging" when several sources 
are present as in the middle and lower panels of Fig. 5. 

The middle panel assumes particle emission from two sources (QT and IS). It is useful 
to compare this plot with those of Fig. 3. One can see that, indeed, particle measured 
in the Ar + Sn reactions can be attributed to the contribution of two main sources, 
a QT and an IS. 

Finally, the bottom panel assumes particle emission from all three sources (QT, QP 
and IS) as parameterized in [18]. Comparing this plot with the Ar + Al data in Fig. 
2 one can notice that: 

1) The QP limiting region is much more defined in Fig. 5 as compared to the ex- 
perimental data (Fig. 2). This is due to the inclusion of more forward angles in Fig. 
5. 

2) The lines in Fig. 5 converge much earlier to a common value of Psource, as Ecm 
is increased, while the experimental data (Fig. 2) present a much larger spreading in 
the intermediate Ecm region. Thus, a broad range of impact parameters leads to a 
"continuum" of sources, while the standard three-moving sources parameterization 
performs an averaging over several sources. 



^ Low and high energy cutoffs in the Psource values are defined by the possibility to extract 
simultaneously Ei, E2 from the invariant spectra for the chosen pair of angles. 
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Particle Energy, MeV 

Fig. 1. Proton, deuteron and triton Lab energy spectra produced in 61 MeV/nucleon ^^Ar 
+ ^'''Al collisions, for some angles as listed in each panel. The data are the open and filled 
circles. The upper panel represents protons (points) and Firestreak model calculations (his- 
tograms), normalized to data at 78°. Coalescence model calculations (Sec. 5) for deuterons 
(middle panel) and tritons (bottom panel) are plotted as histograms. The (PN/d^dE scale 
is presented in ij.b/{sr ■ MeV). 
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Fig. 2. The "effective" source velocity Psource versus the kinetic energy (in the source frame) 
of the protons (full symbols), deuterons (big open symbols) and tritons (small open symbols) 
emitted in 61 MeV/nucleon ^^Ar + ^^Al collisions for 4 pairs of angles. The solid lines are 
drawn to guide the eye. 
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Fig. 3. The "effective" source velocity Psource versus the kinetic energy (in the source frame) 
of the protons (full symbols), deuterons (big open symbols) and tritons (small open symbols) 
emitted in 61 MeV/nucleon ^^Ar + ^^^Sn (upper panel) and ^^Ar + ^^^Sn (lower panel) 
collisions for 4 pairs of angles. The solid lines are drawn to guide the eye. 
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Fig. 4. Summed over impact parameter (with bmax= 12-7 fm) the yield function, Y{ri), the 
velocity, /?, and internal energy per particle, t, are plotted against the projectile fraction 77 
for 61 MeV/nucleon ^^Ar + ^'''Al collisions. 
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Fig. 5. The "effective" source velocity Psource versus the proton kinetic energy in the source 
frame, as reconstructed from the energy spectra assuming particle emission from single 
(upper panel), two (middle panel) and three (bottom panel) moving sources. The parameters 
of the Maxwell-Boltzmann distribution in the source frame are listed in the figure (QP is 
quasi-projectile source, QT is quasi-target source and QF is the intermediate velocity source, 
in the notation of Ref. [18]). 



20 



